Surface states in zigzag and armchair graphene nanoribbons 
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This paper presents electronic spectra of zigzag and armchair graphene nanoribbons calculated 
within the tight-binding model for 7r-electrons. Zigzag and armchair nanoribbons of different edge 
geometries are considered, with surface perturbation taken into account. The properties of surface 
states are discussed on the basis of their classification into Tamm states and Shockley states. In 
armchair nanoribbons surface states are shown to close the energy gap at the Dirac point for certain 
edge geometries. 
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I. INTRODUCTION 

Graphite-based materials are the subject of inten- 
sive research for both fundamental and practical rea- 
sons. Fullereneai, carbon nanotubea^ and graphene 
nanoribbonsii^ have unusual electronic properties, which 
allow the design and fabrication of nanoelectronic sys- 
tems with parameters unavailable to conventional elec- 
tronics in this size range-. 

Graphene, a single layer of graphite, is composed of 
carbon atoms arranged in a two-dimensional honeycomb 
lattice. Covalent a bonds between hybridized sp2 or- 
bitals form the skeleton of a graphene sheet. Electronic 
transport in graphene is, in principle, based on elec- 
trons from pz orbitals forming delocaHzed tt bonds. The 
tight-binding approximation (TBA) commonly used for 
a rough description of transport properties of graphite- 
based materials^J^ only takes into account the nearest- 
neighbor hopping of TT electrons. 

The unit cell of the honeycomb lattice comprises two 
lattice sites. Consequently, two energy bands of graphene 
are obtained in the tight-binding model with hopping 
limited to nearest neighbors. The bands touch at six K 
points in the Brillouin zone. Exactly between the bands, 
the Fermi energy passes through the K points. The lin- 
ear character of the dispersion relation near the K points 
results in electrons behaving as massless fermiona^ de- 
scribed by the Dirac equationiS; hence the K points are 
referred to as Dirac points. As a two-dimensional system 
with a specific topology of the band structure, graphene is 
characterized by the occurrence of non-zero Berry phases 
of the electronic wave functionii. This brings about both 
the integer and the fractional quantum Hall effect, pre- 
dicted theoreticallyi^*12. and verified experimentall y^"^'^^ . 

Carbon nanotubes and graphene nanoribbons are 
quasi-lD systems formed by rolling or confining, respec- 
tively, a graphene sheet along the direction defined by 
the so-called chiral vector Vc- The assumed boundary 
conditions imply the quantization of the wave vector 
component kc describing the propagation of an electron 
wave in graphene in the direction of Vc- The corre- 
sponding dispersion relation i?„(A:||) for the propagation 



along the nanotube/nanoribbon comprises a number of 
branches resulting from the quantization of k^. This 
implies quantized conductance in carbon nanotubes and 
graphene nanoribbons. At low temperatures electronic 
transport in quasi-lD graphene structures presents spin- 
related effects, such as the Coulomb blockad o^^i^^ or the 
Kondo eSec^F^, due to the electron-electron interaction. 

In carbon nanotubes all the dispersion branches He 
within the projection of the 2D graphene dispersion re- 
lation on the direction of kc, as implied by the peri- 
odic boundary conditions and the lack of surface^. In 
graphene nanoribbons, due to their finite size in the di- 
rection of Vc, the quantized wave vector component kc 
can take on complex values. The corresponding states 
localize at the nanoribbon surface (edge)22i>2ii2^ and their 
dispersion branches i?„(/cj[) He beyond the projection of 
the graphene dispersion relation on the direction of kc- 

The orientation of kc in the Brillouin zone also deter- 
mines the position of the Dirac points in the dispersion 
relation En{k\\). For kc oriented along the T — M direc- 
tion Dirac points occur at fcy = ±|7r . Rotating kc to 
the T — K direction results in the Dirac points shifting 
to fcji = 0. 

The transport properties are determined by the elec- 
tronic structure near the Fermi energy. As in graphene 
the Fermi energy passes through the Dirac points, the 
orientation of kc with respect to the Dirac points, and 
the circumference/length of the nanotube/nanoribbon 
are decisive for its metalHc or semiconducting charac- 
ter. By suitable adjustment of these parameters kc can 
be quantized in a manner which implies some dispersion 
branches passing through the Dirac points. This results 
in the energy gap being closed by branches of bulk states. 

In nanoribbons the energy gap can be closed by an- 
other mechanism as well. Metallic character of the sys- 
tems (by assuming simple TBA model) can be due to 
the occurrence of surface states^i between bands of bulk 
states. In zigzag graphene nanoribbons (with kc oriented 
in the T — M direction) strongly localized surface states 
occur near the Fermi energy. The role of these states 
is especiaHy important in thin nanoribbons, with a high 
density of surface states in a relatively wide gap between 
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the upper-band bulk states and the lower-band ones. No 
surface states have been shown to occur in graphene 
nanoribbons of pure armchair edge structure (fee along 
the T — K direction) in the tight-binding approach^°. For 
kc orientations between the T — K and T — M directions, 
dimers of carbon atoms (characteristic of the armchair 
orientation) as well as single atoms (typical of the zigzag 
orientation) occur at the nanoribbon edgesSA. Surface 
states have been shown to localize at surface atoms typi- 
cal of the zigzag edge and their density to have a peak at 
the Fermi energy. Interestingly, the nearly fiat dispersion 
branches £'„(fc||) in the vicinity of i5 = and the con- 
sequent high density of states imply spin-polarized edges 
when Hubbard repulsion is taken into accoun t^^'^^ . Spin- 
polarization results in half-metallic properties of zigzag 
nanoribbon^. It is worth to note that the result of the 
electron-electron interaction is the opening a small gap 
at Fermi level both for zigzag and 'metallic' armchair 
nanoribonsSS.. This processes were not taken into account 
in the presented paper. 

Two factors determine the occurrence of surface states 
in the system considered in the model approach: bro- 
ken translational symmetry and surface perturbation. 
These two factors provide the basis for the distinction, 
used by some authors2^i21i2^^ of two categories of sur- 
face states, referred to as Shockley2£ states and Tamir^ 
states. Shockley states are surface states predicted to oc- 
cur at a non-reconstructed surface conserving the chem- 
ical composition of the bulk. The conditions of existence 
of Shockley states at a non-perturbed surface can be for- 
mulated on the basis of the symmetry of the system. The 
surface perturbation, when taken into account, tends to 
impair the localization of surface states of this type. In 
contrast to Shockley states, Tamm states necessitate a 
surface perturbation to occur. Real surfaces are always 
perturbed due to reconstruction and chemisorption pro- 
cesses. However, in spite of the concurrent occurrence of 
broken translational symmetry and surface perturbation, 
the discussed classification can be of use for studying the 
conditions of existence of surface states in the system. 

Many studies of graphene nanoribbons in the tight- 
binding approach assume no surface reconstruction and 
no significant effect of hydrogen passivation of the surface 
on the energy of surface carbon atoms. The main purpose 
of the passivation, performed in experimental studies, is 
to saturate the sp2 bonds of surface carbon atoms. As a 
result, zigzag edges are composed of alternately arranged 
surface and bulk carbon atoms (having two and three 
neighbors, respectively), while armchair edges consist of 
alternately placed dimers of surface atoms and dimers of 
bulk atoms. 

A free radical in the form of an extra hydrogen atom 
or a methyl group can be added or at the nanoribbon 
surface. This addition results in modified hybridiza- 
tion of surface carbon atom orbitals, which become sps 
hybrid a^^i^^ . There is also possible to consider saturation 
of some danghng carbon bonds by methylene groups^^ 
(not by hydrogen). This delivers additional 7r-electrons 



coming from each —CH2 group. Processes in which the 
chemical composition of the surface is modified can be 
considered, too. It can by done by substituting boron 
atoms or nitrogen atoms in place of carbon atoms-- or 
by passivation of surface by another radicals: —F, — 0, 
—OH—. These effects can be taken into account in the 
tight-binding model by modifying the energy of a sur- 
face carbon atom32. and its neighbor hooping, or, more 
radically, by excluding a carbon atom from the surface 
(or adding an extra one) . This will modify the geometry 
of the nanoribbon skeleton in which electrons propagate. 
In the resulting model each surface atom will have only 
neighbor with one bulk carbon atom. 

In this paper surface states localized at graphene 
nanoribbon edges are determined in the tight-binding 
approximation with nearest-neighbor hopping. Different 
model geometries of armchair and zigzag nanoribbons are 
considered, with surface carbon atoms having two bulk 
neighbors or a single bulk neighbor. The surface pertur- 
bation is modeled by a modification of the surface atom 
energy. The effect of energy gap closing by bulk states or 
surface states in armchair and zigzag nanoribbons and its 
dependence on the edge geometry are discussed as well. 

The ensuing part of this paper is organized as follows: 
Section II discusses the considered graphene nanoribbon 
structures. Section III presents the method employed in 
the calculations, the results of which are presented and 
discussed in Section IV. The study is summed up in the 
closing Section V. 



II. MODEL STRUCTURES 

The structures depicted in Figs[T] and [2] represent the 
skeleton of graphene nanoribbons with 7r-electron delo- 
caHzation. Figures [T] and [2] provide a scheme for the 
description of the electronic properties of the system in 
the tight-binding model; circles and fines represent lattice 
sites and hoppings between them, respectively. 

Our interest is limited to structures without surface 
reconstruction. The geometry of a reconstructed surface 
can differ substantially from that of the nanoribbon bulk 
and practically represent a different phase^^ . Quite chal- 
lenging, a systematic description of such systems can be 
the subject of further studies. 

The nanoribbons shown in Figs. [Hand [2] are cut from 
a graphene sheet in two particular directions, resulting in 
two types of edge termination, referred to as zigzag and 
armchair. In graphene nanoribbons with other orienta- 
tions of the chiral vector the edge structure has many 
steps and terraces^L^^ and can be regarded as a mixture 
of the zigzag structure and the armchair one. 

Let us consider two model edge configurations of the 
nanoribbon skeleton, labeled / and // in Figs. [T]and[2l 
In the type / configuration, each surface (edge) site has 
two bulk neighbors, while in the type // configuration 
each surface site has only one bulk neighbor. The type / 
configuration, with surface carbon atoms sp2 -hybridized 
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Figure 1: Structure of 7r-electron skeleton of zigzag graphene 
nanoribbon. Circles and lines represent lattice sites and hop- 
pings, respectively; large circles represent surface (edge) sites. 
Each surface site has two neighbors (in edge configuration /) 
or a single neighbor (in edge configuration //). Structures 
AI and AH have a center of symmetry, while structures SI 
and 5*// have a symmetry axis. Structures AI — II {AI com- 
bined with AH) and SI — II {SI combined with SII) have 
edges of both types. Dashed lines are limits of rectangular 
unit cells of graphene sheet (periodic in x and y directions); 
gray area represents a cell of graphene nanoribbon (periodic 
in y direction). The nanoribbon cell comprises two slices of 
lattice sites. Surface perturbation is introduced by modifying 
the energy (potential) of surface sites. 



as a result of hydrogen saturation of the broken bond, is 
regarded in the Hterature as that of a non-reconstructed 
pure surface of graphene nanoribbon. However, a free 
radical added to a hydrogen-passivated surface atom 
makes it spa-hybridized and excluded from the skeleton 
of sites participating in transport of 7r-electrons. Regard- 
less of the physical mechanism, discussed in the previous 
section, we can presume that surface sites with a single 
bulk neighbor should be taken into account in the tight- 
binding model. Configuration // represents a Hmiting 
case. In real systems only a part of surface atoms have a 
single neighbor, due to reduced radical-carbon bond dis- 
sociation energy or steric restrictions, which increase with 
coverage of the surface by radicals/citesurf-H-bound-dys. 

Depending of the nanoribbon width, nanoribbons with 
both edges of the same configuration (/ or //) can have 
a symmetry axis (structures SI and SII) or a center 
of symmetry (structures AI and AH). Systems with 
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Figure 2: Structure of vr-electron skeleton of armchair 
graphene nanoribbon. Symbols and classification rules are 
as for the zigzag structure (Fig. [T|). Note the nanoribbon cell 
(periodic in x direction) comprises four slices of lattice sites. 



'mixed' edges (structures AI — II and SI — II) are easily 
seen to have neither a symmetry axis along the nanorib- 
bon nor a center of symmetry. 

A hexagonal unit cell comprising two atoms can be 
defined in an infinite graphene sheet. Repeated unit 
translations of the unit cell atoms produce two sublat- 
tices, distinguished by empty and filled circles in Figs. [T] 
and [21 Note that a different lattice cell is more conve- 
nient in the case of graphene nanoribbon. A rectangu- 
lar cell comprising four atoms reproduces the structure 
of the nanoribbon when translated along and perpen- 
dicularly to its axis. The area of this nanoribbon cell 
is twice as large as that of the hexagonal unit cell of 
graphene. Consequently, the corresponding rectangular 
Brillouin zone (with vertices in M points) is formed by 
folding the hexagonal Brillouin zone (with vertices at K 
points), see Fig. [3l The area covered by the bulk disper- 
sion branches En{k\\) of an armchair or zigzag nanorib- 
bon is determined by a simple projection of the 2D dis- 
persion relation of graphene on the T — K direction or the 
T — M direction, respectively, from the areas delimited by 
elongated rectangular frames in Fig. [31 This operation 
is equivalent to a projection of the folded Brillouin zone. 

In our discussion of the geometry of graphene nanorib- 
bons we have not taken into account surface perturba- 
tion so far. We have only allowed for different orienta- 
tions (armchair or zigzag) of the nanoribbon against the 
graphene lattice and different positions of edges in the 
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III. METHOD 




zig-zag 



arm-chair 



Figure 3: Two-dimensional dispersion relation of graphene 
for upper/lower band; dark areas correspond to energy values 
close to zero. The hexagon with vertices at K points delimits 
the first Brillouin zone. Inside, the rectangle with vertices at 
M points represents the folded Brillouin zone corresponding 
to the rectangular unit cell of the original lattice. The elon- 
gated rectangles are areas of projection of dispersion relation 
on r — if and F — M directions. Bottom and right, graphene 
lattice against edges of armchair and zigzag nanoribbons, re- 
spectively. 



unit cell (configurations I and //). Graphene nanorib- 
bons are often regarded as unrolled carbon nanotubes. It 
is worthy of notice that among the systems under consid- 
eration only zigzag structures SI and AH and armchair 
structure AI can be rolled to form a carbon nanotube. 
These structures have an integer number of rectangular 
unit cells (delimited by dashed line in Figs. [1] and [2]) 
across the nanoribbon. Below we shall attempt to deter- 
mine which of the model structures under consideration 
permit the existence of Shockley states, or surface states 
which do not require surface perturbation. 



In the considered model the surface perturbation will 
be introduced by modifying the energy (potential) of sur- 
face atoms in the 7r-electron skeleton. In Figs. [Hand [2] 
this is visualized by enlarged circles representing surface 
sites. We are going to determine the effect of surface 
perturbation defined in this manner on the existence of 
Shockley states, as well as the conditions in which Tamm 
states appear. 



The method presented below allows to determine wave 
vectors of modes of known energy, propagating in the 
nanoribbon along its axis. First developed for the deter- 
mination of band structure of lattices of mesoscopic quan- 
tum dots2£, this technique was subsequently adapted to 
the calculation of graphene nanoribbon spectra^. 

The following symbols and units will be used below for 
clarity of presentation: The wave-vector components 
and ky will be expressed in units of ^ and respec- 
tively, where 3a and \/3a are dimensions of the rectan- 
gular unit cell of the graphene lattice and a is the bond 
length. The obtained magnitudes of the dimensionless 
components kx and ky in the folded Brillouin zone range 
from — TT to TT. The electron energy X is expressed in 
units of the hopping integral t. The assumed zero energy 
level is the energy (potential) in lattice sites e 



X 



E 



t 



(1) 



Considering the translational symmetry along the 
nanoribbon, let us define a cell comprising M = 2 or 
M = 4 sHces of lattice sites in the zigzag structure and 
in the armchair structure, respectively (see Figs. [T]and 
[21). The tight-binding Hamiltonian can be expressed in 
the form of a block tridiagonal supermatrix in which the 
diagonal submatrices (blocks) represent Hamiltonians of 
isolated cells: 



i.i — 
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Ui.i-l Hi U, 
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Hm I 



(2) 



and the superdiagonal = U and subdiagonal 

Hi+i,i = blocks define the intercell hopping: 



U 



( 







Ui.m\ 







(3) 



) 



Element Hi is the Hamiltonian matrix of the z-th isolated 
sHce of lattice sites, and t/i,i±i is the hopping matrix 
for adjacent slices. Note matrices Ui^i±i are generally 
rectangular, due to the possibly different number Nj of 
lattice sites in different slices. In the adopted units of 
energy all the diagonal elements Hi'"' are zero, with the 
exception of H\'^ and H^'''^' , which take on values: 



Z = 



t 



(4) 



if the extreme sites in the slice i are surface sites (e^ 
denotes energy of surface site. The nonzero elements of 
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Hi and Ui,i±i indicate intersite hopping within the sHce 
or between adjacent sHces, respectively. 

Let us divide the state space into two subspaces, one 
referring to a reference cell and the other to the rest of 
the system: 



(5) 



The Hamiltonian of the whole system can be formally 
written as: 



H = H 



cell 



H° 



(6) 



where U defines hopping between the reference cell and 
the rest of the system, defined by the Hamiltonian 
All matrices are written in a whole space of l^"). 

By using the standard definition of a Green's function: 

Qcell ^ _ jjcellyl 

combined with ([5]) and ([6]) and Schodinger equation: 

H\^)=X\^), (8) 
the following relation can be derived: 

\-^,ell) = G^'"U \-^out) ■ (9) 

C^*^" can be determined by the recursive Green's function 
technique^. 

The relation ^ allows to interrelate the functions at 
the interfaces between the reference cell (Z-th) and the 
rest of the system^: 



I r^cell Tj ,x,/— 1 



,1+1 
1 I 



/-icell J J 



(10) 



where columns of the matrix represent wave functions 
for the z-th slice in the Z-th cell. The system of equations 
(flOl) can be put in the matrix form: 



where 



Ti 
Ti 

T2 



^i+i 



= T2 
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By applying Bloch's theorem: 



(+1 



(12) 



(13) 



we obtain a generalized eigenvalue problem, which al- 
lows the determination of wave vectors fcy of nanoribbon 
modes of known energy X: 



Ti 



i-i 

AI 



(14) 



The solutions of the eigenvalue problem l|14p include 
evanescent modes, characterized by complex values of fcy . 
As the system is infinite along the nanoribbon axis, these 
solutions must be discarded as non-physical. Note the 
eigensolver finds solutions corresponding to modes prop- 
agating (or evanescent) in each of the opposite directions, 
and thus with opposite signs of fcy . 

To determine the localization of states in the direc- 
tion perpendicular to the nanoribbon axis we must check 
whether the assumed state energy X for the determined 
value of /cy is within the energy gap of the projection of 
the graphene dispersion relation. The following condi- 
tions of existence of surface states are obtained for the 
zigzag orientation: 



\X\ > l + 2cos(ifcj,) 



(15) 



for states below the lower band and above the upper one, 
and: 

|A:| < 1 - 2cos(ifcj^) for > |7r, 

|A:| < 2cos(ifcj,) - 1 for \ky\<lTr (16) 

between the bands. For the armchair orientation the con- 
ditions become: 



|X| > ^5 + 4cos(ifc,) (17) 
below the lower band and above the upper one, and: 



\X\ < -^3 - 2sin(i7r - |/c,) - 4sin(ifc, + iyr) 

for kx > 0, 



\X\ < •>/3-2sin(i7r+|fc^)-|-4sin(ifc^- iTr) 



for fcx < (18) 



between the bands. 



IV. RESULTS 

Before proceeding to the presentation of results, let us 
discuss the effect of the nanoribbon width on the num- 
ber of dispersion branches and the number of solutions 



Table L Number A'^ of dispersion branches in spectrum and 
number M of modes at fixed energy value in nanoribbons of 
different width n, orientation and edge geometry. 
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Figure 4: Energy spectra of zigzag graphene nanoribbons of 
width n — 10. Plots obtained for geometries 5*/, SI — // 
and SII are shown in respective columns. Spectra in suc- 
cessive rows correspond to growing absolute value of surface 
perturbation: Z = 0, —0.5, —1.5 and —3. Thick lines, com- 
posed of calculation points, represent dispersion branches cor- 
responding to each mode; spacing between points in flattened 
branches is due to flnite energy step in the calculation pro- 
cedure. Grey area represents projection of the 2D graphene 
dispersion relation on the direction of k^^ (cf. Fig. [Sj. Dots 
in surface state dispersion branches correspond to ky and X 
values assumed in the wave function plots in Fig. [U La- 
bels: Is and 2s mark one or two branches of Shockley states 
respectively. 
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Figure 5: Profiles of surface state wave functions plotted 
across a zigzag nanoribbon for the first slices of sites in the 
nanoribbon cell (cf. Fig. [TJ. The states are indicated in 
dispersion spectra plotted in Fig. |4] by black dots. 



of the generalized eigenvalue problem. As the measure 
of nanoribbon width let us assume the maximum num- 
ber n of sites in a slice, allowing for possible slice shift 
(e.g. for structure All in Fig. [l]we assume n = 5). The 
number of dispersion branches and the number of solu- 
tions of the eigenvalue problem for different nanoribbon 
widths are specified in Table [J The number of disper- 
sion branches is easily seen to be equal to that of sites 
in the nanoribbon cell (gray area in Figs. [l]and[2])- The 
number of solutions of the general eigenvalue problem de- 
termines the maximum number of modes propagating in 
the nanoribbon at a fixed value of energy X: the number 
of propagating modes is the total number of solutions of 
the eigenvalue problem minus the number of non-physical 
solutions, corresponding to evanescent modes. 

Figure H] shows spectra obtained for zigzag nanorib- 
bons. Only spectra of structures SI, SI — // and Sll 
are depicted, those of structures Al, Al — 11 and All 
presenting no significant differences. The spectra were 
calculated for nanoribbon width n = 10, which is re- 
flected in the number of dispersion branches: 20, 19 and 
18 in structures SI, SI— 11 and Sll, respectively. Shown 
in rows corresponding to growing surface perturbation, 
the spectra prove symmetric with respect to the sign of 
surface perturbation: its reversal, Z —Z, results in a 
spectrum reflected with respect to X — 0. 

The gray area in Fig. [4] represents the projection of 
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Figure 6: Dispersion spectra of surface states in zigzag 
nanoribbons of width n = 50 and edge geometries SI, SI — II 
and SII with surface perturbation Z = —0.5 and —3. 



the graphene dispersion relation on the direction of k^. 
Branches beyond this area correspond to surface states. 
Shockley surface states are seen to occur at both types 
(/ and //) of non-perturbed surface (Z — 0). For the 
type I edge configuration Shockley states appear in the 
band gap for \ky\ < |7r. If the gap is wide, the value 
of the imaginary component of the wave vector is high, 
which impHes strongly localized states. This reduces the 
interaction between opposite edges of the nanoribbon and 
causes a gradual degeneracy of even and odd states in the 
symmetric structures SI and SII (cf. Fig. [5]). For ky — 
±7r Shockley states at the type / surface are locaHzed at 
a surface site and isolated from the rest of the system, 
hence their energy X — Z. 

As the surface perturbation increases, Shockley states 
shift towards the gap Hmit and progressively lose their 
localization. At the same time, Tamm states emerge 
from the bands. In the case of strong surface pertur- 
bation Tamm states are seen to occur below the lower 
band, above the upper one, and between the bands, for 



\ky 



< 3'^ or \ky\ 



> in structures SI and SII, respec- 



tively. Extreme surface perturbation values cause surface 
atoms to become isolated. This is equivalent to edge con- 
figuration / becoming configuration //, and vice versa. 
This evolution is seen in the dispersion relations obtained 
for Z = and for Z — ~3: with the two branches of 
strongly localized surface states discarded, the spectrum 
of structure SI resembles that of structure SII. The 
mechanism of evolution of the SII spectrum into the 
SI spectrum with growing perturbation is identical. In 
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Figure 7: Energy spectra of armchair graphene nanoribbons 
of width n = 5 and geometries SI, SI~II and SII (columns), 
plotted for increasing absolute value of surface perturbation 
Z = Q, -0.5, -1.5 and -4 (rows). 



structure SI — II discarding the surface sites is equiv- 
alent to exchanging the edge structures, which leads to 
the same spectrum. 

Figure [5] shows wave function amplitudes for the first 
sHce of lattice sites in the cell. The assumed values of 
energy X and wave vector component ky are indicated 
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Figure 8: Dispersion spectra of surface states in armchair 
nanoribbons of width n — \h and geometries SI , SI — II 
and SII with surface perturbation Z — —0.5 and —4. Dots 
in surface state dispersion branches correspond to ky and X 
values assumed in the wave function plots shown in Fig. [9l 



Figure 9: Profiles of surface state wave functions plotted 
across an armchair nanoribbon for the first slice of sites in 
the nanoribbon cell (cf. Fig. [2|). The corresponding states 
are indicated in dispersion spectra plotted in Fig. [8] by black 
dots. 



by dots in surface state branches in Fig. [H The wave 
function profiles testify the locaUzation of each state at 
one surface (in the asymmetric structure SI — II) or at 
both surfaces (in the symmetric structures SI and SII). 
States in the symmetric structures SI and SII are easily 
found to represent even and odd combinations of states 
localized at a single surface of type / or // in the SI — II 
structure. The only important features refiecting differ- 
ence betvifeen S and A structures is the symmetry of wave 
functions. However, the conclusion concerning locaUza- 
tion on the surface of type / or // are the same in both 
cases. In order to prove a localization on both surfaces 
for structures AI and All one have to analyse wave func- 
tions profiles on a series successive slices. Each of them 
is locaHsed on one edge only but the direction of locali- 
sation changes alternately. 

Figure [6] presents surface states in nanoribbons of 
width n = 50 (bulk states are hidden for clarity). In 
a nanoribbon of this width the surfaces are practically 
isolated from each other, which implies a virtual degen- 
eracy of surface states even near the Dirac points, where 
the gap is narrow. The surface state spectra of systems 
of geometries SI and SII correspond to those of a semi- 
infinite graphene sheet of edge structure / or //, respec- 
tively. The spectrum of structure SI ~ II is a superpo- 
sition of those obtained for SI and SII. 



Results of a similar investigation for armchair nanorib- 
bons are depicted in Fig. [7l The presented spectra re- 
fer to nanoribbons of edge geometries SI, SI — II and 
SII and width n — 5. In each of the studied structures 
the spectrum comprises 18 dispersion branches. Each 
surface brings two surface atoms to the nanoribbon cell 
(gray area in Fig. [2|), which implies that at most four 
surface state dispersion branches can occur in the energy 
gaps. Worthy of notice, and only characteristic of arm- 
chair nanoribbons, are energy ranges and 
— \/5 < X < — 1 in which the gap is closed at any value 
of kx. This means no states localized at the nanoribbon 
surface can occur in these energy ranges. 

The spectrum obtained for structure SI indicates that 
no surface states occur at the type / surface without 
surface perturbation {Z — 0). However, Shockley states 
are seen to occur at the type // surface (cf. the spectra 
obtained for structures SII and SI — II at Z — 0). Two 
and four Shockley states occur in structures SI — // and 
SII, respectively, due to the number of type // surfaces 
in each (one in SI — II and two in SII) . 

The spectra of armchair nanoribbons have the same 
symmetry with respect to the sign of surface pertur- 
bation Z as the spectra obtained for zigzag structures. 
Only negative values of Z are assumed in this investi- 
gation. Increasing the absolute value surface perturba- 
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n=4 n=5 n-6 




Figure 10: Dispersion spectra of armchair nanoribbons with 
energy gap closed at Dirac point by width adjustment. The 
surface perturbation is absent {Z = 0). Note the gap is closed 
by bulk states in structures AI and SI, and by surface states 
in the other structures. 



tion causes Shockley states to merge into energy bands, 
with the concurrent induction of Tamm states. When 
the surface perturbation is strong {Z = —4), surface 
sites in each structure {SI, SI — II and SII) become 
isolated, and the system resembles structure SI. This is 
evidenced in the spectrum by the detachment of four sur- 
face state branches from the lower band (note in struc- 
tures SI and SII the detached branches are virtually 
double-degenerate and overlap). 

Figure [8] shows spectra of surface states in armchair 
nanoribbons of width n = 15 (bulk state branches are 
hidden for clarity). Except for the immediate vicinity of 
the Dirac point, the surfaces are seen to be well isolated 
from each other, as evidenced by the close degeneracy of 
surface states in structures SI and SII. By comparing 
the spectrum of structure SI — // with the spectra of SI 
and SII it is easy to determine the type of surface at 
which each surface state localizes in the SI — II system. 
This becomes even more clear if we compare the wave 
function profiles depicted in Fig. [9l The conclusions as 
to the parity of the states with respect to the nanorib- 
bon center in the symmetric structures SI and SII and 
their relation to the wave functions in the asymmetric 
structure SI ~ II are the same as in the case of zigzag 
nanoribbons. 

An interesting effect, which we would like to empha- 
size, is the energy gap closing around the Dirac point. In 
AI nanoribbons the gap is closed by bulk states when the 
nanoribbon has width n = 3to -|- 1, where m = 1, 2, 3, . . . 



Table II: Nanoribbon width n corresponding to energy gap 
closing at Dirac point in armchair nanoribbons of different 
edge geometry; m is an integer. 



structure 


AI 


AI-II 


All 


SI 


SI-II 


SII 


n 


3m+l 


3m+2 


3m 


3m 


3m+l 


3m+2 


Bulk states can 


also close 


the 


energy 


gap in 


an SI 



nanoribbon of width n — 3m (which means two slices 
in the nanoribbon cell comprise 3m sites and the other 
two 3m — 1 sites). In the other structures the gap is 
closed by surface state branches. The nanoribbon width 
and the number of atoms per nanoribbon cell for which 
the energy gap closes near the Dirac point in each of the 
structures under discussion are specified in Table [TTl Fig- 
ure [10] presents spectra obtained for each structure with 
nanoribbon width adjusted so that the effect of gap clos- 
ing can be observed. The dispersion relations plotted in 
Fig. [ini correspond to zero surface perturbation. 

Note the accidental degeneracy due to the intersec- 
tion of dispersion branches in structures SI and AI (with 
both edges of type /) is partially ehminated in structures 
SII and All (with both edges of type //) and does not 
occur at all in structures SI — II and AI — // (with 
mixed edges). If the nanoribbon is not wide enough, the 
ehmination of the degeneracy due to the intersection of 
dispersion branches will result in the generation of energy 
gaps between bulk dispersion branches. 

V. CONCLUSION 

We have determined the energy spectrum of zigzag and 
armchair graphene nanoribbons in the tight-binding ap- 
proximation for electrons. Two model edge configura- 
tions, / and //, have been considered, with surface sites 
having two neighbors or a single neighbor, respectively. 
A surface perturbation has been allowed for and mod- 
eled by a modification of the surface site energy (poten- 
tial). Shockley states, or surface states which do not 
require surface perturbation, have been found to occur 
at edges of both configurations in zigzag nanoribbons. 
In armchair nanoribbons, only the type // edge config- 
uration permits the occurrence of surface states of this 
category. The generation of Tamm states by the surface 
perturbation and its effect on the occurrence of Shock- 
ley states have been examined as well. For both struc- 
ture, a sufficiently strong perturbation destroys Shockley 
states. Localization of Tamm states is, in general, im- 
proved with increasing surface perturbation. The only 
exceptions are Tamm states in the gap between bands for 
armchair structure. They disappear in the one of bands 
for a large surface perturbation. We have also determined 
the conditions of energy gap closing at the Dirac point 
in armchair nanoribbons. In armchair nanoribbons with 
both edges of type / the gap can be closed by bulk states; 
otherwise (in structures II or I — II) the gap can only 
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be closed by surface state branches. 
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